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introduction 

Optical coherence tomography (OCT) is a relatively new technology, which is capable of 
noninvasive micron-scale resolution imaging in living biological tissues. Recent OCT 
research has focused on developing instrumentation appropriate for imaging in clinical 
settings (e.g. in ophthalmology [1,2], dermatology [3] and gastroenterology [4]), on 
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resolution improvements [1], real-time imaging [5 J, and on functional imaging such as in 
color Doppler OCT [6]. A major impediment to the use of OCT for quantitative 
morphological imaging is image distortions which may occur due to several mechanisms, 
including nonlineariries in the reference or sample scan mechanisms, noiv-telecentric 
(diverging or converging) scan geometries, and the refraction of probe light in the sample. 

There are many applications of OCT for which accurate Iwo-dimensional reconstruction 
of complex features in layered media is required. One application of imaging biometry is in 
the ophthalmic anterior segment [2], in which not only linear distances but also curvilinear 
surfaces, interfaces, and enclosed areas must he accurately obtained. In this application, 
ultrasound biomicroscopy (UBM) is a valuable tool for the diagnosis of appositional angle 
closure, which is a risk factor for progressive trabecular damage, elevated intraocular 
pressure and acute angle-closure glaucoma [7]. Since OCT is non-contact, imaging the 
anterior chamber (AC) angle with OCT [2] greatly improves patient comfort, and may 
significantly reduce examination time, thereby enabling screening of populations at risk. An 
additional advantage is the substantial resolution improvement from -50 to <15 urn, which is 
achieved with OCT. 

Current-generation real-time OCT systems employ depth-priority scanning with the 
axial scan implemented using a rapid-scan optical delay (RSOD) in the reference arm [5,8]. 
The rapid axial scan is readily implemented using resonant scanners, however the resulting 
sinusoidally varying delay axiaiiy distorts the resulting OCT images. By contrast, the lateral 
scan is typically slow (4-32 Hz) and may be implemented using closed-loop galvanomclric 
optical scanners without causing image distortions. The use of non-telecentric scan patterns 
is often necessitated by non-planar sample configurations, e.g. imaging the convex surface of 
the cornea or the concave surface of a hollow organ or tract. Non-contact imaging, one of the 
primary advantages of OCT, also leads to significant image distortions due to retraction of 
the probe beam at the interface between air and smooth surfaces such as the cornea or liquid 
accumulations in Internal organs. Image distortions due to retraction may also occur at 
internal smooth tissue index boundaries such as the cornea-aqueous interface in the eye. This 
paper describes algorithms for numerical correction of the three primary unavoidable image 
distortions in real-time OCT: nonlinear axial scanning, non-telecentric lateral scanning, and 
retraction at smooth boundaries. The first two of these algorithms arc demonstrated to be 
corrected in real-time prior to image display, while the third requires interactive user 
definition of smooth tissue boundaries and was thus implemented in post-processing. 

Methods 

Image transformations 

In principle, image transformations can be implemented using cither forward or backward 
mapping approaches. In the raw (uncorrected) OCT image, we define X* and y' to denote the 
coordinates across and along A-scans (single depth scans), x and / are the corresponding 
coordinates in the target (corrected) image (Fig. 1). Pixel positions arc abbreviated as 
P'= (*'/) and P= (x,y). A forward mapping approach (P= /(P)) calculates the target 
position for a given raw image data point. Since the target position P will most likely faQ 
between target pixels, algorithms must be applied to distribute the value of the raw image 
pixel onto several neighboring pixels. Since the density of this distribution is in general non- 
homogeneous, the brightness of the target image must be renormalized, a significant 
computational expense. A backward mapping approach (P J - F(P)) avoids this disadvantage 
by mapping each target pixel to a location in the raw image, then using simple interpolations 
between surrounding pixels in the raw image to obtain Che target pixel value. Furthermore, 
the number of transformations is at the optimum of one per target pixel, another advantage 
compared with forward mapping. If (he same backward transformation is applied to all 
images, it can be implemented with lookup table to achieve real-time imaging [4,9]. 
Transformations can be cascaded, using the output of one transformation as the input for the 
next. 



£932 - 515.00 US Received Feb/ua/y 25,2002; Revised April 23.2002 

(C) 2002 OSA 6 May 2002 / VqL 10, No. 9 / OPTICS EXPRESS 398 



PAGE 24/41 1 RCVD AT 3/2912005 8:54:31 AM [Eastern Standard rime] f SVRill^TO-ffXRF-l/0 1 DNIS:87M306 * CSfD:216 621 6165 1 DURATION (mn>ss):13-22 



03/29/05 09:02 FAX 216 621 6 165 



RENNER OTTO 20 



-> USPTO GENERAL ®025 




Fig. L A) Correcriort of nonlinear axial scanning B) Coordinate system of ihe data, acquired 
(raw image). C) Target coOrdmaio system hi a homogeneous medium; D) Rcftacrive 
hrtcrfeces in the sanjple alter ray paths. Abbreviations: At full amplitude of deplh scan, t\: 
axial dury cycle of acquisition, D: distance cm mos to image of scanning pivot (ISP), w. width, 
tfi depth of acquired image, primed coordinates in raw data, unprimcd in target image, P Wld 
corresponding points in target and source image, scan angle, <p~**: extreme scan angle 
in the vertical center of the image, L*' distance from P or Pi to the ISP, if , Lf- partial 
dismnccfi in inbomogeneoiu sample with the indices of refraction n, and n^L = U+Li+L* 

Nonlinear axial scanning 

The motion of Ihe resonant scanner can be well described as sinusoidal, simplifying the 
backward transformation. In a standard high-speed OCT system the carrier i3 generated by 
Ihe scan mirror motion [5,8], Therefore, the nonlinear motion also limits the usable duty 
cycle T| if no tracking filters arc employed. Within this duty cycle, the scan depth d is 
acquired. With this knowledge, the forward transformation for the axial scan coordinate y 
can be written as (Fig. 1 A) 



(1) 



The total scan depth A can be calculated from 



(2) 



(3) 



Therefore, the backward transformation can be written as 



(4) 
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since the horizontal scan is linear. The rransforniations to correct non-telecentric scanning 
and refraction can use the target image of this transformation as the raw image or (he 
transformations cun be analytically cascaded. 

Non-telecentric scanning 

To correct for geometric image distortions due to non-telecentric scanning, the coordinate 
systems for the raw and target images are defined in terms of the scan geometry (Fig. IB, C). 
Wc assume the OCT image is formed using a sample arm probe Lhai incorporates a rotating 
mirror scanner [2]. The origin for both systems is in the center of the image, x' is linear with 
acquisition time and thereby with the lateral scanning mirror movement. The center of the 
field of view (FOV, liaving width w and depth d) is a distance O from the scanning pivoi 
(SP) or an image of the sample pivot (ISP), depending upon the internal details of the sample 
arm probe. For positive or negative D, scans diverge or converge, respectively. For 
tclccenlric scans D approaches infinity, but for any real system, this can always be 
approximated by a large D. 

The target pixel position P can also be defined in polar coordinates (fl L)> with the scanning 
angle 9 in the FOV and the distance L from ISP to P (fig. IB, C). For a homogeneous 
sample without refraction (denoted by subscript /?), paod L are given by 

tp„ (x, y ) » arctan(x l(D - y )) , ( 6 ) 

and 



^(x.y)=£>-V* z *(D-y) 2 . <7> 

The scanning angle to reach the extreme of the FOV at the center plane is given by 
tPma* =q> h (w /2.0). In the rectangular array of acquired data, the scan angle qt in the FOV is 
linear with position x? <p'(x\y) = 2x" qt^ fw = xVO, while the distance between ISP and 
P is L , (x* 1 /) = 0-/- Since p= & L = L\ and 9^ = > the complete backward 
transformations are given by: 

Refraction correction 

In order to correct for refraction, the refractive indices of sample layers must be known and 
the interfaces between them must be identified (Fig. ID). Then, fundamental principles can 
be applied to correct for beam deflection and to transform optical into physical pathlcngths. 
In the following, we assume k layers of different media in the image area, having indices of 
refraction n, to with the probe beam entering the sample from the top of the image in air 
{no - 1). The interfaces /*(x) between layers of differing index arc assumed to be smooth 
functions. 

A forward transformation for refraction correction could use Snell's law to calculate the 
target pixel for each raw data pixel, by propagating the incident beam sequentially through 
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the sample layers and keeping track of optical path length until the correct distance is 
reached However, the raw image is distorted by ihe scan geometry and refraction on 
previous interfaces. Therefore, defining the normal on the interface and the incident angle to 
apply SnclTs law becomes very complicated in raw image coordinates. For the backward 
transformation, Ihe interfaces can be defined distortion-free in target coordinates, avoiding 
this difficulty. However, for the latter transformation, Snell's law cannot be applied since the 
route of the beam through the sample is not known a priori. A solution can be obtained, 
however, by applying Fennat's principle, which slates that light rays will propagate such that 
the optical path between source and target locations is minimized. Assuming the sample 
beam passes thowzh the points P, at the interfaces of ihe layers to reach P (Fig- 1 D), the total 
optical pathlenglh is the sum of the distance from ISP to Pi (Eq. (7)) and Ihe optical path 
between subsequent Pfs to Pi 

l(p„...p„p)= u W+&»,|»v^|+n.M- (10) 

By varying the location of P, along the interfaces, a mmimum of L can be found, satisfying 
FennaVs principle. Assuming the paths of the different probe beams do not cross within the 
f OV, a unique solution exists for the P\. After iterarively solving for the values of P h the 
complete back transformation can be written as 

^F,(Pi Pa.P) = W)> < U > 

/=F y (P v ...,P k .P)=D-L(P 1 P k ,P). 0.2) 




Fig. 2. Versa tHo lateral scanning hand held probe with * divergent scan used in these 
ex perim ents. Abbreviations: CL: collimaticn ksS, SM: scanning minor, RL: relay lens, OL-- 
Dhjecdve leas, CG: calibration grid with a line spacing of 3 17 pn% translated stepwise oxiaUy. 
f^i focal length of lens XX. 

Experimental Setup 

The high speed OCT system employed in this letter was based on a previously published 
design [5], employing a Fourier-domain rapid-sranning optical delay (RSOD) in the 
reference arm incorporating a resonant scanner oscillating aL 2 kHz. In the sample arm, a 
versatile lateral scanning hand held probe was used, having a divergent scan and a focal 
depth of 1 1 mm [2] (Fig. 2). Two relay lenses imaged the scanning mirror to the aperture of 
the final objective lens. Both relay lenses were achromats to minimize optical aberrations. 
This scanner was designed for small focal spot size, compact design, and sufficient working 
distance to allow for non-contact imaging of the skin or anterior segment of the eye. The 
image size of 3.77 mm width and 4 mm depth was determined by imaging a metal grid with 
317 ^im line spacing (LCI Supplies, Weal Chester, PA), translated axialhy in steps with a 
precision linear stage (Fig. 2 right). A broadband (X=1.32 pm, AX?=68 nm, FWHM), higfr 
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power.(22 mW) semiconductor optical amplifier light source illuminated ihe interferometer. 
For improved SNR, we used an efficient, optical circulator-based interferomcier with 
balanced detection [10]. The axial duty cycle T| of the acquisition was 66.7 %. Wc imaged a 
test object and the temporal anterior chamber angjc in the eye of a healthy volunteer. All 
images were acquired at 8 frames/s to reduce motion artifacts. 

Image processing 

The images were preproccssed online m real-lime to correct the nonlinear axial scanning 
[Eq. (5)] and nou-lelecentric scanning [Eqs. (8) 3 (°)]- The Format's principle algorithm 
derived above was implemented using MatLab 5.2 and applied offline to the acquired images 
in the following steps: Since the geometrical distortion due to scan divergence was already 
removed using the algorithm described above, the first interface was distortion free. 
Therefore, the first refraction interface could be semi-automaticalry defined by user input of 
3 10 4 points on the interface and refined by active contours [11], a popular segmentation 
algorithm. Distortion due to refraction a: the first interface was then corrected using the 
Format's principle algorithm. After the correction of retraction at the first interface, the 
second interface was distortion-free and could be defined in the same manner as the first 
boundary. The Fermat's principle algorithm was then applied again for refraction at the 
second interface. 

The calculation of the Fennat's principle algorithm in its basic form requires a numerical 
iteration tor each target pixel- The computation time was greatly reduced by dividing the 
target image into tiles. In this case, the iteration was only calculated for the comer points of 
the tiles. For all the target points within the tiles, a bilinear interpolation was employed to 
obtain the intermediate transformations. Therefore, for tile sizes sufficiently large, the 
effective computational expense to calculate the transformation for a target pixel was 
reduced by a factor of -100. Since this is an approximation, small errors may occur, 
especially if a discontinuity (a strong refractive index change) crosses the tile. We chose a 
tile size of 10x10 pixel, limiting the error to less then 0.8 pixel. 

Results 
Calibration 

First, we confirmed die precise distortion correction of the resonant scanner movement. The 
calibration grid was translated axialry in steps of 100 um, and the position of the central grid 
wire was determined by the intensity centroid. The maximum position error decreased 
through the transformation from 176 um to 4.6 um (Fig. 3 A). 

Next, the calibration grid was used to validate the accuracy of the geometric image 
transformation for the hand-held probe. An overlay of images acquired at different axial 




SMndWBWrpm 



rig. 3. A) D iscrvticn coirccoan of nonlinear axial seaming- B) Overlay of calibration 
Images acquired at diffcrttnl axial positions, eaeh image averaged over 50 frames, 
corrected far nunlinear axial scanning. C) Rcsidn?) error vector* in air (thin lias), whfccm 
correction (dark Hoc), and wirh aberration correction (length x 10). 
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positions (already corrected for nonlinear axial scanning) of the grid is shown in Fig. 3B. The 
steps between axial positions were 500 unu except the top and bottom steps which were 400 
and 300 um, respectively. The region of best focus differed in the lateral direction, being 
closer to the top on the sides of the JrOV. After calculating the cenlroids of the calibration 
points, we used the best fit to obtain the size of the FOV and the distance D between the 
image of the scanning pivot and the center of the FOV, The FOV had a width w of 3.77 mm 
and a depth d of 4.0 mm. D was estimated as 1 L4 mm, slightly larger than the focal length of 
the objective lens. Fig. 3C displays the residual error for each grid point as a vector between 
the true and estimated location. The vectors for the corrected image are magnified ten times 
for visualization. The maximum error vector length in the FOV decreased from 304 urn in 
the raw data to 17 urn following geometric transformation. 

Refraction phantom 

As a sample test object fox correction of multiple distortions, wb imaged a drop of Intralipid 0 
(rv=1.33) on a glass microscope cover slip (ffjpl.51, measured 10 be 977(10) um thick). Fig. 
4A shows three obvious distortions: (1) The boundaries of the flat cover slip appear convex, 
due to the geometric distortion of the diverging scanner. (2) The upper and lower surfaces of 
the cover slip under the drop appear concave due to refraction by the drop. (3) The cover slip 
appears thicker than its physical thickness due to refraction by the glass. The image was 
acquired slightly displaced from the apex of the droplet to avoid saturation of the detection 
electronics by the large specular reflection which occurs there. Following image correction 
with the algorithms described above, both cover slip surfaces were rendered flat with a 
maximum deviation from a line of 10 and 13 um respectively, and the thickness of the cover 
sup was measured from the image as 978(4) urn (Fig. 4B). Due to the highly positive 
curvature at the edges of the drop, two wedge-shaped areas were visible below, where the 
data was spread out by the algorithm. This effect may be relevant in ophthalmic imaging of 
patients with contact lenses. 

Anterior chamber angle of the eye 

After validating the distortion correction techniques using the test object, we applied the 
techniques to correct images of the anterior chamber angle of the eye of a healthy volunteer. 
The width of the anterior chamber angle and the configuration of the peripheral iris, as 
represented by the area between the inner comeo-scleral wall and the iris, are important in 
the diagnosis of primary angle closure glaucoma [7], Visible in Fig. 5, the appearance of this 
angle changes substantially with processing. The correction of the nonlinear axial scan (Fig. 
5B) expanded the center and compressed top and bottom of the image, while the non- 
telccentric correction compressed the top and expanded the bottom laterally (Fig. 5C). Since 
the corneal surface was significantly tilted, refraction played an important role (Fig. 5D). 



The focus of development in OCT research is shifting from simply acquiring an image to 
developing images that can be used for diagnosis. As shown here, image processing can 
significantly reduce tlie amount of geometric errors within OCT images, allowing for 



Discussion 




tf g, 4. A) Uncorrected image of an [nonlipid 0 drop on a covers) ip, B) corrected image 
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Fig. 5. OCT i rouges °t the temporal unterior chamber angle of a human eye, imaged m 
vivo v 8 fpa. A) raw hnnge, B) corrected Cor nonlinear axial scanning (duty cycle 
6G.7 Q additionally corrected for divergent scan geometry and D) final image, 
corrected for refraction at the oir-COrttt interface (/l^ - 1 .36) and at the 
endoihclium-equeous boundary (/7 BgwBUS = 1 .33) 



accurate measurements of distances, angles, and areas. The main sources of distortion 
requiring correction include nonlinear axial . scanning, non-tclecentric scanning and 
refraction. 

In the general case, refraction is a Ihree-di rnensional problem. Since conventional OCT 
acquires images in two dimensions, accurate correction of refraction can only be achieved if 
no significant beam deflection occurs in the laird dimension. In imaging the anterior 
segment, positioning of the B-scan normal to the corneal surface can be readily achieved by 
monitoring the large specular reflection, which occurs at the apex. In the future, two- 
dimensional lateral scanning could be used to allow three-dimensional refraction correction. 

"While we have already implemented correction of geometric distortion in real-time, 
refraction correction was calculated offline. This was due to both the need for interactive 
designation of refractive surface contours, as wdl as computation time in our preliminary 
implementation. Our MatLab implementation allowed the correction of the images within 7- 
8 sec on a 1 GHz Pentium III computer. Implementing the algorithm in C++ will reduce 
processing time by at least an order of magnitude. In combination with automatic interface 
detection, this should allow online correction during the acquisition. 

Spherical aberration can also influence image distortions. Wc demonstrated with exact 
ray tracing simulations that such aberrations can substantially deviate the ISP when higher 
scan angles or different probe optics arc employed than those used for the data in Fig. 2. 
Furthermore, under this condition the scan angle in the image is no longer linear with the 
angle of the scanner. The aberrations in our system proved to be small enough to be 
neglected, but in cases where the aberrations are significant, these distortions can also be 
corrected using a backward transfbmiation- 

The refraction correction described here is highly sensitive to position noise of the 
detected interface, because the slope varies more with noise than the position. This- 
complicates application of this technique to rough surfaces and may limit its use to smooth 
interfaces. 

Conclusion 

In conclusion, we have shown that major sources of image distortions in OCT, 
refraction, scanning geometries and non linearities, can be corrected from several hundred 
micrometers down to 10-20 urn, reaching the diffraction and coherence length limits. With 
this accuracy, OCT has great potential for micrometer-scale biennctry as a fast, noninvasive 
and quantitative imaging device. 
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